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Abstract 

Sum-over-histories quantization of particle-like theory in curved space is discussed. It is 
reviewed that the propagator satisfies the Schrodinger equation respective wave equation with a 
Laplace-like operator. The exact dependence of the operator on the choice of measure is shown. 

Next, modifications needed for a manifold with a boundary are introduced, and the exact 
form of the equation for the propagator is derived. It is shown that the Laplace-like opera- 
tor contains some distributional terms localized on the boundary. These terms induce proper 
boundary conditions for the propagator. This choice of boundary conditions is explained as a 
consequence of a measurement of particles on the boundary. 

The interaction with sources inside of the domain and sources on the boundary is also 
discussed. 



Introduction 



The main goal of this work is to investigate a quantization of a relativistic particle using the 
sum-over-histories approach. We have the following motivations to do such a study. 

The usual approach how to quantize a relativistic particle is the scalar field theory. Yet, it is a 
quantization of a completely different system — a quantization of a continuous field on spacetime. 
It is true that we can identify some states of such system as particle states — states with some 
properties of particles. But is there any other way to make a connection to the particle theory? 
Is it possible to quantize a classical relativistic particle, and does it give predictions equivalent to 
predictions of quantum scalar field theory? 

There exists a candidate for the direct quantization of a particle theory — quantization using 
the sum-over-histories approach. The classical explanation of this approach for usual nonrelativistic 
physics can be found in Q and more technically in [^-^J. A nice non-technical overview for a 
relativistic theory can be found in ||. Beside these classical introductions, this approach has received 
considerable attention in recent years (see for example Q). The new development has led to a 
generalization of this method called generalized quantum mechanics (see [^|-^)). 

In this approach the transition amplitudes associated with the chosen criteria are computed 
by summing over amplitudes of all possible histories which meet the criteria. It is known that 
some of these amplitudes computed for a relativistic particle lead in the special cases (e.g., in flat 
spacetime) to quantities which can also be obtained from scalar field theory. The goal of our work 
is to investigate this correspondence in more detail. 

There are more reasons for studying relativistic particle theory. One of the attempts to under- 
stand the quantization of the gravitational field coupled to matter is to reduce the full gravitation 
theory to a system with a finite number of degrees of freedom and try to quantize this simplified 
system. These reduced theories are called minisuperspace models. It is well known that such re- 
duced system is essentially equivalent to a particle theory in a Lorentzian space with (usually) a 
complicated potential. A common method for the quantization of minisuperspace models is the 
sum-over-histories approach. 

Key feature of our investigation is that we study the particle theory on a bounded domain 
of spacetime, and that we pay an attention to the exact form of boundary conditions. The usual 
approach is a bit generous on this question — the theory is usually formulated on the whole spacetime 
with not always clearly formulated special behavior at infinities. In the flat spacetime such an 
approach is justifiable because there exists a preferred behavior at infinities, but in a general curved 
spacetime we have to be more careful. The question of boundary conditions is usually completely 
ignored in definitions of the path integral. We try to formulate the theory in a more careful way 
and identify its boundary-condition dependence. 

The plan of our work is the following. In the first part of the paper we review sum-over-histories 
approach to the quantization of the particle in curved spacetime without boundary. Equations 
for key amplitudes (the propagator and Feynmann Green function) are derived and their exact 
dependence on the definition of the path integral is shown. 

It the second part we investigate the theory on a bounded domain. On the bounded domain 
we have to modify our definition of the path integral. The new definition leads to a modification 
of the equations for physical amplitudes. Some new distributional terms localized on the boundary 
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appear in these equation and we find that they specify the exact form of boundary conditions for 
the amplitudes. It is shown that the exact form of the boundary conditions depends on the details 
of the definition of the path integral. 

Finally we discuss shortly a physical meaning of the boundary conditions. We argue that bound- 
ary conditions can serve as a phenomenological description of apparatuses measuring particles on 
the boundary of the domain. 

We will continue in this discussion in the following paper where a relationship of different 
boundary conditions will be investigated and a connection with a particle definition in the quantum 
scalar field theory will be shown. 

Sections 1 and 2 contain the main line of arguments, details of computations can be found in the 
appendix ^[ The appendix [b| contains a general overview of the geodesic theory necessary for the 
computations, including not common theory for reflected geodesies. 

Let us note that our study can accommodate a wider range of theories. The theory will be 
parameterized by two signature factors. One of them (the factor n) characterizes a signature of a 
target space metric (the metric in which a particle lives) , and another (the factor v) describes whether 
the theory is physical or Euclidian. For real n and imaginary v we obtain non-relativistic particle in 
curved space, for n and v both imaginary we obtain relativistic particle in curved spacetime and for 
real signature factors we get mathematically better behaving, but non-physical Euclidian version of 
the theory. Usually most quantities are well defined in the Euclidian sector and the definition for 
physical signature are obtained by an analytical continuation in the signature factors. 
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1 Particle in a curved space without boundary 



Space of histories, action and amplitudes 

In this section we shortly formulate the sum-over-histories approach to a quantization of a particle- 
like theory. A similar calculation have been done for example in JTo|,|TT[|. We present our derivation 
here because we will generalize our line of reasoning in then next section to the case when a particle 
is moving in a domain with a boundary. 

As usual in the sum-over-histories approach, the theory is characterized by a space of histories 
and an action. An elementary history in our case is a trajectory — an imbedding of a 1-dimensional 
manifold N (called the inner space) to a d-dimensional spacetime manifold M — and an inner space 
metric h on the inner manifold N. M is equipped with a spacetime metric g and scalar potential y. 

In the Euclidian version of the theory h is positive definite; in the physical version it is negative 
definite. In the Euclidian version we allow the spacetime metric to be also positive definiteEl. 

The whole theory is invariant under diffeomorphisms of the inner manifold N. As usual (e.g. 0]), 
we factorize over this symmetry. If we fix a coordinate 7/ : N — > (0, 1) on the inner manifold N we can 
characterize a class of equivalent histories using a pair [x, r], where a; is a map x : (0, 1) — > M 
and r is a total inner time or a total inner length of N measured using the inner metric h. 

The Euclidian actiono in these variables has the form 



(04) 



+ vtV{x) )dr] ■ (1.1) 



Here v is a constant signature factor distinguishing Euclidian and physical versions of the theory. 
In the former case v — 2, in the later it is multiplied by i. 

In the sum-over-histories approach to quantum theory we can define an amplitude A(H) for any 
set of histories H by "summing" over amplitudes of elementary histories in the set. The quantum 
amplitude is not directly a physical measurable quantity. We need an additional notion of distin- 
guishable or decoherent histories to give a probabilistic interpretation to the square of amplitudes. 
We expect that this notion has the same symmetry as the action and a measure on histories. This 
means that we will be always interested in amplitudes of sets of histories which are invariant under 
the action of the diffeomorphism group. For such sets we can factorize the path integral and elimi- 
nate the reference to the diffeomorphism (e.g. |l2|). In the factorized integral we are summing only 
over variables [x.t]: 

A(E)= J Tlred(r,x) exp(-J(T,x)) , (1.2) 

[i,i]£H 

with a reduced, renormalized measure 9Jt re( j- 
Propagator 

It is useful to compute an amplitude —K(t, X{ |x;) - called the propagator or heat kernel - for the 
set of histories restricted only by positions of end points of the trajectory Xf and X{ in the spacetime 
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M and by fixing an inner time to a particular value r: 

^K{T,X{\xi)= J m F {T,xi\ Xi )[x} exp(-I(r,x)) , (1.3) 

where 7"(xf IxiLis a set of trajectories x : (0,1) — > M with x(l) — Xf and x(0) = Xi and n is a 
constant factorS governing signature of the spacetime metric g. 

Because the set of histories [r, T(xf \xi)] = {r} x T(xf\xi) is a lower dimensional subset of the 
space of all histories, K(r, Xf|xi) is essentially an amplitude "density" on the space R + x M x M of 
values [r, Xf, xj. Therefore we have to expect that the restriction OJI^t, Xf |xi) of the measure 9Jl r ed 
to the space [r, T(xf |x;)], which we call the Feynman measure, depends on r and end points Xf and 
Xi; maybe only in a "trivial" way. 

In other words, an amplitude density of an elementary history on the space [r, T(x{ \x{)] is 

A{t,x) =M F (T,Xf\xi)[x} exp(-I(r,x)) . (1.4) 

It is well known Sfl[Tc|,|ll|| that with the right choice of the measure 9JI f (t, Xf|xi) the propagator 
satisfies the equationQ 

-K(t) = ^F.K(t) , ^(OHS" 1 , (1.5) 

where F is b wave operator fixed by the action and the measure (see later), Q = q^5 is a delta 
distributional on M normalized to the metric volume element g 2 . In other words, K is the exponential 
of F 

X(r)=cxp(-yF).g- 1 . (1.6) 

We have not specified the "right choice" of the measure yet. It can be a very problematic 
task from the pure mathematical point of view. Instead of trying to develop a measure theory on 
infinite dimensional spaces for oscillatory integrals (where the main problem lies), we take the usual 
approach of formal manipulations, and we define the measure by its decomposition properties and 



approximation for small time intervals. The former is given in equation ( |l . 1 1 ) , and the latter is 



given in equation ( 1.20 ) 



The idea of the proof of the relations (1.5) is in proving key properties of the exponential, 



K(r f ) • Q • K( Ti ) = Kfa + n) , (1-7) 

g,K( T ).g = g -^-t + o{t 2 ) , (1.8) 



where T = g • F is the quadratic form of the differential operator F. 
Composition law 



The first condition (1.7) is a composition law for the amplitude K. This law reflects the possibility 
to decompose a history [t,x] into histories [Ti,:Ei] during an initial amount of inner time t\ and 
[n, Xf] during a final amount of inner time n ■ We say that a history [t, x] = [t{, x{\ © [ri, x{\ is given 
by joining of histories [71, Xf] and [t\, xj if 

r = r f + r i , Xf(0)=Xi(l) , x f = s(^±^) , a* = x(^) . (1.9) 

T T 

The actions is additive with respect of joining histories. 



5 



We have a natural decomposition of the set of histories [r, T(x{ \xi)] which defines the propagator 
K(t, Xf\xi) to disjoint sets [t{, T(x{ \x q )] x [T u T(x \xi)] 

[T,T(xf\xi)] = [j [T U T( Xi \x )} x [n,T(x \xi)] . (1.10) 

If the measures on these sets are related by 

Tt F {r,Xf\xi)[x} =m F (T i ,x f \x )[x{] nQ^(x )Tt F (n,x \xi)[xi} , (1.11) 



we get (L7) by a straigthforward calculation. 

The condition (1.11) represents a reasonable assumption of the locality of the measure Wl F . 



Together with the additivity of the action it reflects the rule of the sum-over-histories approach 
to quantum mechanics — that the amplitude of independent (here consequent) events is given by 
multiplication of individual amplitudes. This condition is the first part of our definition of the 
measure. Now we know how to construct the measure 9JI f (t) for some time r from measures 
for shorter time intervals. To conclude the definition of the measure, we need to specify it for an 
infinitesimally short inner time interval. This moves us to an investigation of the short time behavior 
of the heat kernel. 



Short time amplitude 



Now we turn to prove equation (l.S). It can be found in the literature (e.g. |lfj,|llj]), but we present 
it here to show how the measure is actually determined and how the operator T depends on this 
choice. 

We ignore technical difficulties in the definition of the path integral, and we assume that this 
integral has most of the properties of a usual integral in a finite-dimensional manifold. This allows 
us to find the short time behavior for the propagator. 

First we write an expansion of the action for small r 

I(T,x) = -I-l(x)+I (x)+Th(x) + ... . (1.12) 
T 

For the action we are using it means 

I-i(x) = ±- J x a x f3 g af3 (x)dr ] , (1.13) 

•76(0,1) 

l (x)=0 , (1.14) 
h(x) = ^ J V(x)dr) . (1.15) 

rj6(0,l) 

We assume that the measure is slowly changing in r compared to the leading term in the action. 



The dominant contribution to the integral (1.3) comes from an extremum x(x{\xi) of the leading 



term 7_i in the exponent. But the extremum of the functional ( 1.13| ) is clearly a geodesic of the 
metric g. We expand all expressions around this extremum 

x = x( X f \xi) + \frx , (1-16) 

where a; is a tangent vector to the space of trajectories T{xi\x{) at the extremum x(xf\x{). We 
actually need to specify what the addition in the last equation means. It will be done more carefully 
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in a similar situation in appendix [A| (see eq. (A. 5)). But now we are interested more in a qualitative 
answer, so we skip these details here. The expanded integral (1.3) has the structure 



— K (r, Xf \x{) — exp^— l(x{xf \x{))j J WI f *(t, xf\xi) exp ^— ^x ■ 5 2 I-i(x(xf ■ x 

(l + V^ (x odd - terms) + r (:r cvcn - terms) + . . 



(1.17) 



Here 9JI f *(t, Xf \x{) is a leading term in the r and ^-expansion of the measure VJl F (t, Xf\xi) after 
change of variables x — > x. DJl F « is a constant measure on the tangent vector space to the space 
of trajectories T(xf\xi). The actual dependence on x is hidden in higher terms of the a;-expansion. 
As a leading term in the r-expansion, 5DT F * depends on r in a trivial way — it is proportional to a 
power of r. Of course, this statement is formal — the exponent of r in 9Jt F * is of the order of the 
dimension of the tangent space, which is infinite. 

"ce-terms" in the last equation represents terms resulting from the expansion of the action and the 
measure; x odd or £c GVOn suggest that x occurs in these terms in odd or even power. For convenience 
we combined the term rl\ into the prefactor despite the fact that it could be included among terms 
proportional to r. 

The value vl~*\ {x(x{ \xi)) is a well-known quantity called the world function, or half the squared 
geodesic distancaa, 



a(x{ \x{) = vl-i(x(xf \xi)) = i / x(xf\xi) ■ g(x(xf\xi)) ■ x(xf\xi) drj 



We also use the notation 



V(x f \xi) 



<o,i) 



-h(x(x{\xi)) = J V(x(xi \xi)) dr\ 
<o,i) 



(1.18) 



(1.19) 



The integral (1.17) is a simple Gaussian integration. (In fact, one approach to defining infinite- 
dimensional integrals is through the definition of a "Gaussian" measure which in our case would be 
971^* exp(— i x ■ <5 2 7_i • xf). The integration can be performed, at least formally. If the measure 
WI f *(t) has the already mentioned r-dependence, the result can be written 



K(t, Xi\xi) 



(2itvt) 2 



|-A(x f |a;i) \ ao(xf\xi) - olx{x{\x{) + 0{r 2 )j x 



c\p [ -— a(x{ \xi) - ^-V(x[\xi 



where ao(xf\xi) should satisfy 



ao(x\x) = 1 



(1.20) 



(1.21) 



Here A(xf\xi) is Van Vleck-Morette determinant (see ( B.13 )). 

The terms proportional to y/r disappeared during integration, thanks to the odd power of x. The 
par ticu lar behaviour a coincidence limit of the coefficient ao will be needed for a proper normalization 
in ( |l.5| ). To obtain this behaviour we need the mentioned r-dependence of the measure, which can 
be expressed by condition 



m F 4r,x\x) Det 



6 2 I-!{x(x\x)) 

2llVT 



1 



(1.22) 
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I.e., the meassure must satisfy this condition to conlude proof of equation (1.8). 

Coefficients in front of powers of r could be expressed in terms of variations of the action and 
the measure. But because we did not specify the measure precisely yet, we can do it now by fixing 
these coefficients. I.e., we can define the measure 0Jl F by choosing functions ao(x{ \xi) and ai(xf |xi). 

In the following we will prove that a form of the operator T in ( |1.8| ) depends only on the 
coincidence limits of ot\ and the first two derivatives of qp. So we can ignore terms with hi gher 
power of r in cq. ( 1.2C| ). As discussed before, equations (1.8) together with composition law (1.7) 



determines the propagator K, i.e. also all important information hidden in the measure 9Jl F . This 
means that a knowledge of the mentioned coincidence limits concludes our definition of the measure 
and path integral itself. 

Let us note that this argument is some kind of justification of the usual time-discretization of the 



path integral and of a priori choice of the short time amplitude in the form ( 1.20 ). But in principle 
it would be possible to define the measure VJl F in some more compact way and compute exactly the 
form of the functions czq and ot\ in terms of variation of the action and the measure. 

Short time behavior of the heat kernel 



Now we continue with the proof of equation (1.8). We show that for small r the amplitude ( 1.20 ) 



has the desired behavior in a distributional sense. Most of the technical work is done in appendix 
|a| where it is shown that for small r the following expansion holds (equation ( A.S)) 



J*(x)a*(z) A(x\z) exp( <r(x\z))<p(x)ip(z) = 

(2twt 2 J V vt J 

K 1 x.zeM (1.23) 

= ip • Q • i[) — T—ip»C»tp + 0{t 2 ) , 

where tp and ip are smooth test functions and C is the Laplace operator quadratic formlfl 

<p»C»ip = 1 05 (dtp) ■ g- 1 ■ (dVO = 

, (1-24) 

^(v 2 ^) = - / 0*v(vV) . 

M M 

Let us remember that now we are discussing the case of a manifold without boundary, and therefore 
we do not have to worry about boundary conditions for the Laplace operator and integration by 
parts. 

Using this result it is easy to show that 



M 



ip»g • k(t) »q •ip = 

n 



(2ttvt) 2 



r V - <p . ((d^Hdiuao]) . g + g . {[d^a^g 1 ^ d„)) . ^ + 
0(r 2 ) , 



\i(x)ai(z)A{x\z) (a (x\z) - t V - a x {x\z) + C(t 2 )) 



I l/<j- _ 

ex P ( -— a ( x \ z ) ~ -yV{x\z) ) tp(x)tp(z) 
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where we used [V] = V and [oto] = 1. Here the L4] denotes a coincidence limit of a bitensor A, d r A 
and diA are derivatives with respect of the right and left arguments and the bi-distributions d and 
d are derivatives acting to the left and to the rightS. 
If the condition 



[d r a ] = [diao] = 
is satisfied, we see that the propagator K(t) has really the form ([O]) with 

T=C+V , 



(1.25) 



(1.26) 
(1.27) 



I.e., T is a Laplace operator with a potential term which include the original potential V from the 
action and additional parts depending on the choice of the measure. 

A common choice for ao is a power of the Van Vleck-Morette determinant 



a = A~P , 



(1.28) 



which satisfies the condition (1.25). It leads to an additional part in the potential, 



(1.29) 



where R is a scalar curvature of the metric g. 



The condition (1.25) is actually a consequence of (1.21) and an assumption of the symmetry of 



O' 



a (x\z) = a {z\x) 



(1.30) 



It is a natural assumption in the case when the theory is symmetric under trajectory reversal. 
However this condition does not have to be satisfied if there is a preferred path direction as for 
example in the case of interaction with an electromagnetic field. But we will not discuss such a 



situation, and in the following we will assume that the conditions (1.3C) and (1.25) are satisfied 



In summary, we have seen that for small r the propagator K{j) has the behavior given by (1.8). 
If the measure is defined using the decomposition property fll.7| ) and the short time amplitude ( 1.2C ), 
the operator T is fixed by knowledge of the coincidence limits of a\ and the first two derivatives of 
a . 



Feynman Green Function 

For the relativistic particle the inner time is physically undetectable and therefore any physical 
set of histories will include elementary histories with all possible inner times. Therefore we are 
interested in the amplitude called the Feynman Green function —G F (x{ \xi) associated with the set 
of histories restricted only by the initial and final points Xf, x\. We can obtain it from the propagator 
^■K (t, Xf \xi) by summing over all possible inner times r using the measure ^dr 

-G F (x { \x i )= f -K(T,xt\xi)^dr . (1.31) 
n J n 2 

R+ 



Using eq. (1.6) we immediately get that the Feynman Green function is the inverse of the wave 
operator T . 
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Boundary conditions 



In this section we completely ignored the question of boundary conditions using as excuse that we are 
working in a manifold without boundary. Certainly this is correct, if the manifold M is compact. 
But it is also correct in the case of a non-compact manifold with a sufficiently "nice" metric g 
at infinities. In such cases there exists a canonical choice of boundary conditions for differential 
operators used above, and these boundary conditions usually allow us to integrate by parts. But 
problem arise for the relativistic particle when the manifold M is Lorentzian and operators as C, 
T are hyperbolic. In this case the choice of boundary conditions at the temporal infinities plays 
an important physical role, and it is worth further investigation. First let us note that in a special 
situations (e.g. existence of a time-like Killing vector in the distant past and future) a canonical 
choice of boundary conditions still exists. But canonical here essentially means the most natural 
physical choice. In a general spacetime we do not have this special choice, and we have to address 
the question of the boundary conditions. To deal with this problem, in the next section we will 
investigate our theory in a bounded domain Q of the manifold M. 
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2 Particle in a curved space with boundary 



General consideration 

In this section our goal is a better understanding of the physical meaning of boundary conditions 
for the differential operators in the equation for the propagator and Green function. Therefore we 
restrict ourselves to a domain Q in the manifold M which is bounded "in all physically interesting 
directions" . This means that we will investigate the boundary conditions on boundaries which are 
not at infinity. 

If the target space metric is positive definite (i.e. in Euclidian and non-relativistic version of the 
theory) all problems with infinities could be solved by a restriction to a compact domain. But in the 
case of relativistic particle, which we are mostly intersted in, we also allow domains which do not 
have to be compact, we allow the domain SI to be unbounded if we know that its infinity is "safe". 

The situation we have in the mind is the Lorentzian globally hyperbolic manifold with asymp- 
totically flat spatial infinity. In this case we can ignore spatial infinity because it makes a sense to 
restrict ourselves to situations in which spacetime is "empty" sufficiently far in space directions. But 
because of the hyperbolic nature of the evolution equation we cannot ignore boundary conditions in 
the time directions. They represent "initial" and "final" conditions of the system. And we want to 
understand exactly this relationship. 

Therefore in the case of a Lorentzian globally hyperbolic manifold, the typical choice of the 
domain will be a sandwich domain between two Cauchy surfaces. But whole discussion is also valid 
for Euclidian version of the theory restricted to a compact domain. 



Restriction to a domain — naive approach 

Let us start with a straighforward restriction to a domain SI. We want to compute an amplitude 
K q (t, X{ \xi) which corresponds to a set of histories with inner time r, endpoints x\, Xf and which 
wholly belong to the domain SI. We can repeat the derivation of the short time amplitude ( 1.2 0| ) , at 



least for X{, Xi sufficiently far from the boundary, because for small r only trajectories near to the 
geodesic between xt and x\ contribute to the amplitude. 

If we do this calculation, we find that there is a new term in the expansion. As can be seen in eq. 



(A. 7), the smoothed short time amplitude leads to a Gauss integration in a variable Z from a tangent 
space at a point x, and in the case of a space without boundary the integration of odd powers of Z 
disappears. But in the case of a domain with a boundary for a point x near the boundary the Gauss 
integration is not always over the whole tangent space and therefore the integral of odd powers of Z 
does not disappear. As shown in appendix |A| (equation (A. 26)), the correct asymptotic expansion 



of the leading term of the short time amplitude (equivalent of ( 1.23| )) is given by 



(2-kvt) 2 

x,zGM 



1 , * 1 



02 (x)g2 ( z )A(x\z) exp ( a(x\z) ) Lp(x)ip(z) 



VT 



(2.1) 
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where Q[9f2] is a delta bi-distribution localized on the boundary normalized to the boundary volume 
element understood as a distribution on spacetime, 



on 



(2.2) 



and £ is a particular ordering of the Laplace operator given by 



C 



(2.3) 
(2.4) 



Using this result it is easy to show that the expansion of the propagator K Q is 

and T is a Laplace-like quadratic form with potential, 

T T = T = 2 + V . 



(2.5) 

(2.6) 
(2.7) 



The corrected potential V is given again by the expression ( 1.27j ) and we have assumed that the 
condition ( 1.25| ) is satisfied. 

We see that the expansion of the propagator has an additional term localized on the boundary <90 
proportional to yfr. This r-dependence causes a problem because K o (0) is singular on the boundary. 
An origin of the singular term on the boundary lies in our careless approximation of the propagator 
by the short time amplitude (1.20). This approximation is correct only for endpoints sufficiently 
far from the boundary. For points near the boundary we have to investigate the structure of the 
propagator more thoroughly. 



Boundary correction term 



The short time amplitude (1.20) represents the dominant contribution to the heat kernel from 
trajectories near the geodesic joining endpoints Xf and X{. But in the case of a sum over trajectories 
restricted to the domain O there are other dominant terms given by contributions of trajectories 
near extremal paths which reflect on the boundary. 

In general we should take into account trajectories with an arbitrary number of reflections on 
the boundary and compute the dominant contributions from all of them. However, for endpoints 
sufficiently far from the boundary the contributions from the reflected paths are negligible compared 
to the straight geodesic — for small r only short paths contribute to the sum, and any trajectory 



with a reflection on the boundary is too long (see figure 2.1). 

But for endpoints near to the boundary the contributions from the reflected trajectories can be 
comparable with the leading term. For the endpoints near a smooth boundary there exists exactly 
one extreme trajectory x^(x\z)) with one reflection whic h giv es a contribution comparable to the 



contribution from the straight geodesic x(x\z)) (see figure |2.l|) . 



The reflected extreme trajectory a?b(a;|z)) is an extremum of the leading term of the action (1.13) 



with the additional condition that the trajectory reflects on the boundary. Let's denote the point of 



12 



(B(B 



Figure 2.1: Example of extremal trajectories. 

Dominant terms to sum over trajectories are given by trajectories near to extreme tra- 
jectories, possibly reflected from the boundary. If close endpoints are sufficiently far 
from the boundary, the reflected geodesies are longer then the straight geodesic. If the 
endpoints are near the boundary, there is a reflected geodesic with the length comparable 
to the length of straight one. Near the corner there are more reflected geodesies with 
comparable length. 




the reflection b(x\z) and the parameter for which the reflection occurs A r (a;|z) and its complement 
\\(x\z). Clearly the trajectory is a joining of two geodesies 



[r, x h (x\z)] = [X\(x\z)t, x(x\b(x\z))} [X t (x\z)T, x(b(x\z)\z)] 
Using additivity of the action we get the value of its leading term 

/ I \ clef T / _ / I ^ <Tl(x\z) (T T {X\Z) 



Xi(x\z) X r (x\z) 



where, following the convention (B.77), 

&\{x\z) = cr(x\b(x\z)) , a r (x\z) = a(b(x\z)\z) 
The extremum requirement gives us conditions on b and Ai, X\, 

T>a{x\b{x\z)) Da(b(x\z)\z) 



X\(x\z) 



Hi 



X T (x\z) 
A? ' 







(2.8) 

(2.9) 

(2.10) 

(2.11) 
(2.12) 



where D denotes the orthogonal projection of the gradient on the boundary (the gradient with respect 
of the argument on the boundary) . See appendix ^ for more details and other quantities defined on 
the boundary. 

Now we can estimate the contribution from the trajectories near to the reflected geodesic x^{x\z)). 
Using reasoning similar to that used for deriving (1.20), we can write an approximation of the short 
time amplitude associated with the reflected geodesic as 



K h (r, Xf\xi) 



(2itvt)- 



-A* p (x[\xi) P{r,Xi\xi) exp 



-o-bOEflzi) 



(2.13) 



where Ab is Van Vleck-Morcttc determinant associated with the reflected geodesic (see (B.89)). The 
coefficient (3 is an analog of the coefficients «0j ct\i only in this case we have to expect an expansion 
in powers of y/r: 



(3{t,x\z) =0q(x\z) + ^?(3i{x\z)+0{t) 



(2.14) 
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As we will see, the right normalization relative to the leading term K a requires 

(3 Q (x\x) = 1 . 



(2.15) 



We did not bother to write down a potential term, because terms of order 0(t) are negligible in the 
approximation we need, as can be seen in the calculation in appendix |A|. We also already anticipated 
an arbitrary power of the Van Vleck-Morette determinant, similarly to the choice (1.28). 

Fixing this short time amplitude (i.e. specification of coefficients p and f3, or more precisely its 
coincidence limits as we will see below) together with amplitude (1.20) concludes the definition of 
the path integral in the domain with a smooth boundary. 



Short time behavior of the heat kernel 

Next we proceed to derive the short time behavior of the propagator. Again, the technical work is 
done in appendix [X[ where it is shown that for small r we have the expansion (see eq. (A.46)) 

ip»g • K h (r) •Q»ij} = 

.1 



f T-\ — if • Q • ijj - 
n V 2-7T 



(2.16) 



T 2^' 



(ii^fc + /3-terms)Q. V + C(r^) . 



The /3-terms contains coincidence limits of the first two derivatives of the coefficient /3 on the bound- 



ary, and the exact form can be found in (A. 56). k is the trace of the external curvature (B.5C) and 

are defined in ( |2.25|) . The normalization ( |2.15|) ensures that the ^/r-term in Kb cancels 
exactly with such a term in K a . 

So, if we add both dominant terms we get 

K k (r) = K (t) + K h (r) , (2.17) 



(2.18) 



where T^. is the quadratic form of Laplace-like operator with the boundary conditions given by the 
choice of /3 coefficient^] 

(2.19) 
(2.20) 
(2.21) 



— — Qk 

d a .(x9" la/3 e).d 



(XV) 



del 



e k = - 



i+p 



k + /3-terms ) Q 



where x is the characteristic function of the domain f2. Here Tn is the standard quadratic form 
which appears in the action for non-interacting scalar field with a potential evaluated on the domain 
f2 and 0^ is a bi-distribution loca lized on the boundary. 

The short time behavior ( 2.18 ) together with the composition law again prove the heat equation 
for the propagator: 

v 



-g . K k {r) = -T- k . K k {r) (2.22) 

The quadratic form T k identifies what kind of boundary conditions the propagator K k satisfies. 
If we compare T k with the operator T (all derivatives act to the right) we find 

T k =T + dT k , (2.23) 
dT k = dT A -Q k , (2.24) 
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where the bi-distribution dJ-^ is essentially integration of a value of the left argument and a momenta 
of the right argument over boundary of the domain 

ip • <£F g • rp = J if n a d a ilj . (2.25) 

To see clearly what kind of boundary conditions the propagator satisfies, we write the heat equation 
in the following way 

-G»k k (T) = ~lF»Kk(T) + ~dr i »K k (T) . (2.26) 

The solution of the heat equation is smooth for non-zero time r. Therefore the left-hand side is 
smooth, as well as the first term on the right-hand side. The second term is localized on the boundary 
and therefore it has to vanish. So we find 

dF i »K k (r) = . (2.27) 

Or, if we define maps t£ and 7r which assign a value ip = (£• <f> and a momenta ir = n • (f> on the 
boundary to a spacetime function <j>, we can writeo 

dFa=£-7r , (2.28) 
B k = ((^-k + P-teTtas) qh) ■<£ , (2.29) 

and the boundary conditions get the form 

(z + + ^-terrxis)^ • K k (r) = . (2.30) 

We see that K k satisfies Robin-like boundary conditions. 

Let us summarize. We have found that the propagator given by the sum of amplitudes over 
histories in the domain Q, with fixed endpoints and inner time r is a solution of the heat equation 
with specific boundary conditions. The boundary conditions depend on the definition of the path 
integral through the coincidence limits of derivatives of coefficients (3 in the short time amplitude 
( 2.13| ). In general, they are Robin-like conditions with a non-degenerate coefficient in front of 
momentum. 



Green function 

In the case of a relativistic particle the inner time is an unphysical quantity, and all physically 
distinguishable sets of histories should contain histories with all possible inner times. Therefore we 
will compute the amplitude associated with the set of histories with fixed endpoints but without, a 
restriction on the inner time. As before, we will call this amplitude the Feynman Green functionta 

-Gl(x\z)= f -K k (T,x\z)^-dT . (2.31) 
n J n 2 

Using the heat equation and initial conditions for the heat kernel, the integration gives 

T- k • Gf = 5 . (2.32) 

So, is the inverse of J 7 ^ and restricted to smooth sources it satisfies the same boundary conditions 
as the propagator K k . 
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Amplitude of particles emitted by a source 

Let us compute the amplitude zj^ (J) of a set of histories which end at a given point x and are 
emitted by a source described by a spacetime dependent amplitudeQ nJ. We will call it the one- 
particle amplitude. It is clearly given by 



It satisfies the same boundary conditions as the Feynman Green function. 

We will interpret the boundary conditions as a consequence of the fact that we have not allowed 
particles to start on the boundary. More precisely, we have allowed the smooth source to be non-zero 
up to the boundary, but we have not allowed an emission of particles from the boundary comparable 
to an emission from a finite volume. 

We can ask why some particular boundary conditions means that no particles are emitted from 
the boundary. What about different boundary conditions? Why is the choice of the conditions 
above special? We are touching a question of what kind of particles we are dealing with. What 
does it mean that no particles are emitted (or absorbed — for scalar particle the meanings are 
interchangeable if we do not distinguish initial and final parts of the boundary) . 

First we have to realize that the statement "no particles on the boundary" has to be interpreted as 
a result of a measurement on the boundary. We have to arrange apparatuses on the whole boundary 
which are sensitive to particles, and when all these devices measure no particle we can speak about 
no emission or absorption. Clearly this is very complicated global measurement. It depends on an 
exact arrangement of experimental devices on the whole boundary and on an interaction of particles 
with devices. We have hidden this dependence in the definition of the path integral through the 
non-specified /3-terms. Therefore we see that we cannot expect a unique canonical meaning for 
the statement "no particles on the boundary" . Only if we specify the kind of measurement we are 
performing we do have a meaning for this statement. And necessary information about experimental 
devices can be phenomenologically characterized by the choice of boundary conditions of the type 
we encountered above. 

In the following paper || we will discuss what interaction of particles with boundary leads to 
different boundary conditions. 

So, the different boundary conditions correspond to different type of detections of particles on 
the boundary and in the following paper we will see that it can correspond to different definitions 
of vacua in the standard particle theory - in the scalar field theory in curved spacetime. 

Emission from the boundary 

Of course, we can ask what is the amplitude to find a particle at a point x if we allow an emission 
from the boundary. Let us assume that the amplitude of the emission from the boundary is given by 
a density nj on the boundary manifold, which we call the boundary source. The amplitude Z^p (r; j) 
associated with the set of one-particle histories which are emitted by this boundary source and end 
in time r at a point x, can be written using the boundary propagator 



The boundary propagator propagates between points inside of the domain and boundary sources. 
It has the character of a function on the domain 17 in the left argument and the function on the 
boundary manifold dfl in the right argument. 

Clearly, the boundary propagator satisfies a composition law similar to ( |l.7| ) 



(2.33) 



Z^(r;j)=K^r)-j 



(2.34) 



K]?(T)=K k (T-e).g.Kj(e) 



(2.35) 
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We can take a limit e — > and get 

K?(t) =K k {r)*K k , (2.36) 
^^0.^(0) . (2.37) 

We see that the amplitude is given by the propagator K k (r) with no emission from the boundary, 

and by the boundary term K k which "translates" between the space of sources on the boundary 
and amplitudes in the domain. Similarly, if we sum over all possible inner times we get 

zl 1) (j)=G£.K k .j . (2.38) 
The boundary term K k is a zero-time amplitude, so it is straightforward to estimate it. The 



short time amplitude approximation similar to (1.20) for the boundary propagator is 
= 7^7^ / sHx)4,{x)j{y)A{x\y) exp(-— a(x\y))(l + 0(V?)) 

Therefore, for zero inner time we get 



(2.39) 



K k = !& . (2.40) 

It means that the emission from the boundary is equivalent to the emission of particles inside of the 
domain but with a distributional source d J = ]• <£. with support on the boundary 

Allowing both boundary sources and sources inside of the domain, we find that the one-particle 
amplitude is 

4 1} (J,9J) = Gf »{J + dJ) = 4> k {J + dJ) . (2.41) 

A careful discussion of the distributional character of introduced boundary sources, differential 
operators and Green functions shows (see |l3|]) that it satisfies the equation of motion in the expected 
form 

F k • 4>k{J + dJ) = J + dJ (2.42) 
with boundary conditions fixed by the boundary source 

dJ = dF k »4> k {J + dJ) . (2.43) 
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Conclusion 



In this paper we studied sum-over-histories quantization of relativistic particle on a bounded 
domain of the spacetime. We showed that we have to modify a definition of the path integral by 
adding terms corresponding to paths reflected on the boundary of the domain. Such contributions 
can be dominant in the short time approximation near the boundary in addition to the usual 
dominant contributions from the straight geodesic. They compensate other terms localized on 
the boundary which arise from restriction of non-reflected paths to the domain. They also specify 
the exact form of the boundary conditions for the propagator and Green functions. We found 
that boundary conditions have Robin-like form and their exact form depends on the details of the 
definition of the path integral — a non-uniqueness is hidden in the specification of /3-coefficients in 
the short time amplitude (2.13). 

We interpreted the specific boundary conditions as a consequence of an interaction with appa- 
ratuses localized on the boundary which allow us to define a notion of particles. Because of the 
non-uniqueness of the definition of the path integral we do not have a uniqueness in the definition of 
particles. In this correspondence the boundary condition can be also viewed as a phenomenological 
description of the specific kind of particles. 

In the next paper || we will discuss a freedom in boundary conditions in more detail. We will 
show that inclusion of some additional interaction with measurement apparatuses on the boundary 
can lead to general boundary conditions which include the conditions corresponding to a definition 
of particles in the quantum scalar field theory. 
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A Asymptotic expansion of the leading term in 
the heat kernel 



Vector space 

In a vector space V equipped with a positive nondegenerate quadratic form g a simple Gaussian 
integration gives 

— T / S Hx) 2 Hz) lp (X)^(Z)exp(-^-(X^Z)-g-(X-Z) 



(2-KVT 

' x.zev 



f gi(XMX) [ zi{Y)i>{X + ^Y)e W (~Y-g-Y 
(2tw)2 J J v 2^ 

xev Yev 

-±-r I f(i W i)^ T ![3 ai ...^](i) x 



(A.l) 



xev 

X h J 9 h (Y)Y a i...Y°«° exp(-i-y. 5 .Y) = 

rev 



m! 

meN 



Here y>, ^ are test functions, is the constant volume element of the metric g, v is a number 
(Re i ^ 0) and £ m represents a bi-distribution 

tp.£ m .iP = J &h[(9 la0 d a d r^] = J &h[(g- lal3 d a d0r<p\ . (A.2) 

v v 

It can be viewed as an "m-th" power of the Laplace quadratic form C associated with the metric g 
on the vector space V. 

So we can write the asymptotic expansion for small r as 

^— rQ k (X ) Q i(Z)^p(-^-(X-Z).g.(X-Z))= h{- r iT Cm ■ 

Here we allowed the metric g to be Lorentzian — this case can be obtained by analytical continuation 
in the phase factor n which characterizes the signature of the metric. 

In the case of the vector space the expression on the left side of eq. (A. 3) is the heat kernel 
of the operator C Because the right side is formally the exponential we see that the expansion is 
exact. To say more about convergence it is necessary to specify functional spaces on which all the 
operators act and we will not do this here. But see e.g. Q for some details. 
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Manifold without a boundary 

Now we would like to find the expansion of the similar expression in a general manifold M without 
boundary. More precisely, we want to expand 

A(x\z)gi (x)q* (z) exp a{x\z) (A.4) 



{2-kvt) 2 
for small r. 

We smooth both arguments with test functions <p, ip and note that for a small r the integration 
over x and z is dominated by a diagonal x ~ z thanks to a(x\z) « for x s» z. Therefore for a fixed 
x we can restrict integration over z to a normal neighborhood of x. In this neighborhood we can 
change variables z — > Z with 

z = u x {sfrZ) nx + V^Z , (A.5) 



where u x (eZ) is a geodesic with an origin x and initial tangent vector Z (see (B.l)). This is the 
exact meaning of "adding" of a vector to a point as mentioned after equation ( 1.16| ). 

The Jacobian associated with this change of variables is given by Van-Vleck Morette determinant 
(see flR43| )) 

(u x 1 *gi)(V^Z)=rhH^[Z]A- 1 (x\z) , (A.6) 

where (x) [Z] is u nders tood as a constant measure on the target vector space T X AI . After a change 
of variables, using ( B.16| ), expanding i/j and performing a Gaussian integration, we get 

A(x\z)g2 (x)q^ (z)(p(x)ip(z) expf — — a{x\z) 



(2-kvt) 2 

x,z£M 



2»(x)ip(x) / Qz(x)[Z]x 
xeM zet.m 



_ (A.7) 

(E jf^*a,...a t (l) ^...Z afc ) exp (- — Z ■ 9 {x) ■ Z) = 



feeNo 



,^(^ + ^ 5 - la ^ 2 ^ + 0(r 2 )) 



A/ 



Using ( B.27 ) we have -02 = Vd0, so we get 

— ^- T A(x|z)^( a ;)^(z)exp(-- ( T(.x|z)) = ^ - ^£ + 0(r 2 ) . (A.8) 

Half line 

Next we will investigate the simplest case of the manifold with the boundary — half line R + . We 
assume it is equipped with a special coordinate r\ which selects a measure and derivative 

u = dn , M = fid , d=S- ■ (A.9) 

or] 

We define bi-distributions of the m-th derivative 

w d {m) • ip = uj • (d m ip) = J uj{d m ip) , d {m) = d {m)T , (A.10) 
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and 

d = d^ , d = d {1) . (A.ll) 
We can define also a boundary delta bi-distribution V as (the boundary is one point now) 

yV»1p = (<pi[)) | boundary • ( A - 12 ) 

Integration by parts can be expressed by the relation 



d»M + M»d = -V , 

d {m+1) . m + {-\) m M • a (m+1> = - d {m ~ k) • v • a (fc> 



/c— 0,...,m 



Next we define quadratic forms of powers of the Laplace operator 
C = -M»d {2) , 2 {rn) = (-l) m M • d {2m) 

2 = -d {2) . m , 2 {m) = (-i) ro a< 2m > . x 

< £=i(£ + £) , £ > < m > = i(£ <m> + 2 {m} ) 



fe+i=2m 



(A.13) 



(A.14) 



£ = d • M. • d . 

The symplectic form on the boundary is 

dC = 2-2= -d»V + V»d . (A.15) 
It is straightforward to check that 

2<m) _ 2<™> = (_!)"» ^ a (fe> .2?.a <0 . (A.16) 

fe,(eN 

fe+(+l=2m 

Now we prove the following expansion for small r 

\- HTpryH , (A.iT) 

- Z-j m! I 2 J 

m£No 

( ir+1 (^^^..(O 

4^. r(m + |) V 2 / 2 



Here ^ is a complex number such that £ has non- negative real part. Strictly speaking, the following 
derivation needs a positive real part; but for an imaginary value of v the relation can be obtained 
by limiting procedure. Because only a combination the tv appears in the equation, we drop v in 
the following derivation — it can be easily restored by inspection of the r-dependence. 
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Clearly, the second equality follows from integration by parts (A. 13). To prove the first one we 
smooth it with test functions tp and ip and get 



'2-KT 



exp 



(-^-0 2 ) v(Mt) KMC) 



1 






-1- 


1 





(A.18) 



«6<0,e> 



J rfC^(0^(C)ex P (-^(e-C) 2 ) 



{G(e,oo> C6R + 



for some e <G R + . 

For a small r the exponential suppresses any contribution except from £ « £. Therefore for 
small e only small values of £ and £ contribute to the first term of the last equation. We can rescale 
variables by factor yfr and expand <p and ifi at zero and obtain 



= J d£ J dC^(C)exp(-_L(£-C)' 



£e(0,e) CSK+ 
T 



27TT 



/el 



(A.19) 



where we have used the definition of special functions R; (C.l). Properties of these special functions 
are summarized in appendix |c[ Using equation ( C.llf ) the last expression gives 



_L J2 <p< k \0) ^\0) x 



kd£N 



x (-l)'+*R fc+z+2 (0) + 



2tt 



E 



(2m- 1)!! / I 



l\k\ ^ k + l-2m 

me N 
2m<l 



lGmAvf) 



fc+/-2m+l 



e i(^)"' R — (-*) 



=0,...,fe 

£ r^(2^r(^±i + l))"i £ (-1)^^(0)^(0) 



/=0,...,m 



c fe+/+l 



(2m)!! ^ kl U (k + I + 1) 

m6N v ; fc,/eN v ; 

le 2 



+ exp(--i-)0(r) 
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(2m-i-l 



>(0) ^ <0 (0) 



T m 1 
meN v ; !=0,...,2ra-l 

S (-i)',«—>(0)^(0) + 

mGN V ; V '-n 



1=0,. ..,2m 

le 2 



(A.20) 



+ exp(--i-)0(r) 



Here we have used relations ( |C.10 ) and (C.14). 

In the second term of the expression (A. 18) we change variables C ^ V = "^77 (£ — C) using again 
the fact that only the contribution from £ ss £ is not suppressed by the exponential. Next we expand 
tp around r\ = and obtain 



= J d( J dc^(0v>(c)ex P (--^(e-c) 5 



1 

v/27rr 



E (2m)!' 



m=N 



E 



fe(e,oo) 



d^(£)V <2m> (0 



7= E tI _E ^jr J «(|'M-*)^^ 



fcGN /=0,...,fe 



(e,oo) 



(2m)!! 



^(0v> (2m> (0 - E 



m=Nr 



(2m)!! 



(A.21) 



Here we used ( C.15 ), ( C.ll ), ( |C .16 ) and ( C.14 ). For small e we can expand <p and ip^ 2m ^ about zero 
in the second term of the last expression. Performing also an integration by parts in the first term 
transforms the last expression to 



me N 



mGN V ;=0,...,2to-1 

- e e t i„;^' +1 ) ^'>»*""' , "'°» 

/ 1 e 2 \ 
+ exp(-- 7 )0(r) . 



(A.22) 
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Substituting equations ( |A.2C| ) and (A. 22) to equation (A. IS) and ignoring exponentially sup- 
pressed terms exp(— \^t)0(t) we obtain the desired relation ( A. 17 ). 
Next we will prove another expansion for small r 



\I2-kvt + <; / V iVT I 

4, V 2 J rfs±l) n + m + 1 ( n + m ) 



riGN V 2 7 m,fcJeN 



where u> is some smooth function. As a corollary for w — 1 we get 
1 



exp(-J-(e + C) 2 ) M(OA»(0 = 

= y __I n^hw.p.^D^ (A.24) 

fe.ZGNo V 2 T V 



m£N fe,iGN 

fc+2+l=2m 



E E 3 w .i>.3 w 



mGN v 2; fe,ieN 

fc+;=2m 



As before, we drop the factor v in the proof because it can easily be restored from the t de- 
pendence. Again, we smooth the relation with test functions ip and ip. Thanks to the exponential 
suppression, only small values of £ and £ contributes to the integrals. Therefore we can rescale £ 
and C by y/r, expand ip, ip, and u> about zero, and get 



± / ^(C)^(^)exp(^(£ + C ) 2 )^dC = 

= ^£^ + ' £ ^W'>(o)*«(o)> 

1 /• crn+kf-m+l . i 



5.CR+ 

-Z^ v2 J r(^±i) ^ n + m + l ^ W (UJ ' A.25 

k+l+va—n 

what proves the relation ( A. 23 ). Here we have used the integral ( |C.lg| ). 
Manifold with boundary — no reflection contribution 

Now we find an expansion of the contribution to the short time amplitude from the trajectories near 
the geodesic without reflection in the domain Q with boundary dQ. We will prove for small t 



(2itvt) 
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2 (x)j}*(z)A(x\z) expf-— <j(x\z) 



(A.26) 
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As usual, we will be proving a smoothed version of this relation — we multiply the expression 
by test functions <f{x) and ip(z), and integrate over x and z. Thanks to the exponential suppression 
the only non-trivial contribution is from x ~ z. Therefore it is sufficient to prove the relation locally. 
Clearly, for ip and ip with support inside of interior of the domain f2 the boundary does not have any 
influence and the relation reduces to the case without boundary. Therefore we will investigate only 
the case when (p and ip are localized near the boundary. Thanks to locality we can also, without 
losing generality, assume that the test functions are localized on the neighborhood U C O of the 
boundary with topology R x 5f2 on which the geodesies normal to the boundary do not cross. In 
such a neighborhood we can use the method described in appendix [b| and change the integration 
over a neig hborhood to integration over the boundary dfl and geodesic distance from the boundary 
(see ( |B.60| )) 



x— >x,£ , x = w(x.£) , 

y ' (A.27) 

z^2,C i y = w(z,() ■ 

The Jacobian associated with this change of variables is 

fl*(2) = deq*(a,0 = de?(2,0q* ■ ( A - 28 ) 



Here we use the convention ( B.62 ) — we denote the spacetime dependent object A{x) in variables 
x, £ and tensor indices moved to the boundary as Changing variables we get 

gi (x)g^ (z)A(x\z) exp[ — — a(x\z))ip(x)rp(z) = 



(2ttvt)2 J \ VT 



(2ttvt) 2 

2,26<9f2 

£,CeR+ 



<%<K j(x,0^(S)j(z,C)^{z)A(x^\z,C) >< 



x exp 

VT 



V VT I 



n 



q 2 (2) — rr / d£d((p(x,£) x 

{2irvT)2 J 

2em i f f ceR+ , i \ ( A - 29 ) 

x 77; / q 2 (2) A(S|2)?/'(2,C) exp( a(x, f| 2, C) + l(x, £| z, ()) . 

(2itvt)2 J \ VT J 

Here we have defined 

i (s, i% c) - info, A(£ ;51f : - 3 % o) , (a.so) 

V A(x\z) J 

where A(x\z) is the Van-Vleck Morette determinant of the metric q on the boundary manifold. 

Exponential suppression ensures again that the only contributions comes from i»2. So we can 
change variables x,z — > y, Y" 

y = 2 , FeT^J! , 2 = i> g (y7Y) , (A.31) 



with the exponential map v^(Y) defined as in (B.l) but on the boundary manifold. The Jacobian 



for this change of variables is given by an equation similar to (B.43), only with the Van-Vleck 
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Morette determinant A (2, y) denned using the metric q on the boundary manifold. We use covariant 
expansions 



fc=2,3,. 



fc=N 



fc=N 



(A.32) 
(A.33) 
(A.34) 



with coefficients given by expressions ( p!99| - pl0^ ) , ( |B.104| |B.106| ), and ( pT27| ). Expanding the 
exponential, gathering all terms up to order 0(t), and performing a Gaussian integration over Y 
lead to 

2 

dtd( exp(-£-(C-0 2 ) ~mt,0 (^o(y^,C)+ra;i(y;e,C)+0(r 2 )) 



(27wt) 2 



with 



£0 = A(y, e|g, C) J (y, (9) (Det„ a , 2 (y; e, 
= (j(^)A(y,£|y,C)J(S,C)) 5 , 

wi(j/;£,C) = v<p(y>€) $(£/>C) (-g?o,4M«/«A S-'o^ u do;2 X + ^0,2^ S , 2 



(A.35) 

(A.36) 
(A.37) 



1 // 7^ 



1 



r, „ — lct/3~-lK,A T " -1. 

ff ; 3(ira Co,3k.A/3 CT ,2 °0,2 °0,2 — 77 ' 0,1 fx "0,3 kAi/ Cq 2 u 02 



-1 ixv -~_-\ K \ 



(A.38) 



The equality in ( A.36| ) is proved in (B.76). Applying the expansion in (A. 17) in (A.35) and using 
(A. 6) we get 



2 (y) </?(yMy) 



yen 



1 / v 
n V 2-7T 



q 2 (y) v{y)^{y) 



yen 



dvv{y)l{mri) (4^2 i" 1 ^'^) 3 



4n 2 



j" 1 (jw ) r ' r ' + wij (f/; ??, rj) 



(A.39) 



yedn 
»+ 



+ 0(rS) 
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where C) (or / r '(£,C)) means derivative of a function / with respect of the left (or right) 

argument. 

We see that we have proved already the desired expansion up to order y/r. Now we proceed to 
prove it in the order r. We split the integrand of the last term to pieces and compute each of them. 
First we note that (see (B.101)) 



[o-o,2 J = q , 



(A.40) 



where by a coincidence limit of a function depending on two real parameters we mean [/] (77) 

f(r),r)). Using flB.lOTp and flB.108|) we get 



and therefore 



i[(lnif']=lnfc , 



(A.41) 



(A.42) 



Here k is a trace of extrinsic curvature and K is a square of the extrinsic curvature (eq. (B.51)). 
Next, using ( B.27 ), ( B.lOoj ), and ( B.102] ), transforming the connection V to V (eq. ( B.64 )) and 
performing an integration by parts gives 



(A.43) 



Substituting for [<?o.fc] ancl Uo,k] m remaining terms of [toi], a straightforward long calculation gives 



g [<To,4fii/K.Aj q q +-xyiQ,2nu\q 



1 r„ 1 r~ n ~-l n\ If? "IT? 1 "-lww 

1 r. t r„ n ,-1/iv„-1c<)3„-1kA 1 



J- r„ -I r„ -i „-i /xj/„-i ctp„ 

g [0-0,3 Ati/aJ [Oo,3kA/3J q q q 

1 

12 



Oo,3/iKctJ |O0,3i>A/3j q q q — ~['0,l/i| [Oq,3k.Ai/J 9 9 



(A.44) 



= v\ — k' + — K 2 k 2 

1 6 12 8 



Putting together ( A.42 - A.44 ), and integrating by parts we find 



yedn 



v J drj J q 

r+ aa 
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n on 
^ f • (^Cd - -(dCd + dCd)^j • tp = -- f C » ip . 

This concludes the proof of the expansion in ( A. 2(f ) . 

Manifold with boundary — reflection contribution 

Finally we will prove the last expansion we have needed in the main text: 

— £ h p (x\z) (J(t,x\z) exp( — -<r h (x\z)) = 

I 2 \ VT / 



(A.45) 



(27Wr) 2 



t-J^Q-t~ (dTd + dT d ) -r V - (^pfc + /3-terms) + 0(r* ) 



(A.46) 



n V 27r 



2 2 



Similarly to the previous section, we smooth this expression with test functions, perform a change 
of variables (A. 27) and consequently (A. 31), and use the covariant expansions (A. 34), 



^,2^12,0=30(2.^12,0 + ^^(2^12,0+0^) , 
similarly for Pi(x\z); and 

fc=2,3,... 

ibO,cis,v^;o= E ^W^^^O^---^ , 

fc=No 

where 



(A.47) 



(A.48) 
(A.49) 

(A.50) 

A(2|2) 

This leads to a Gaussian integration in the variable Y € T^Sfi in which only leading terms in 
expansions survive and we get 

Q^(x)q?(z) <p(x)ip(z) ^ p (x\z) (3{t,x\z) expf — — a h (x\z) 



■'•;|--o = + ( ^2,0 Abl " P( f,'ff ,C) j(g,0 



(2twt) 



yedn 



(2-kvt)- 



F q%) / d^C ex P (-|-(e + C) 2 ) ^(y,0^(y,C) 



(A.51) 



(8; (A)(5; £,&,0 + v^i (5; + 0(r) 



where 



jb(y; = ^(5; C) j (y, q%) f Det„ a b0 , 2 (y; £, C) 



(J(5,0 A+^^d^C) j(S,0 



A)(9;£,&,0=3oM9,0 > 



(A.52) 
(A.53) 
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with jb and (3q depending analyticaly on their three real arguments. (3i is defined in similar way as 
(3q. Here we anticipate that Ab and f3 can have more complicated analytical dependence on £ and 



C,. The equality in (A. 52) follows from (B.98). Using the expansion (A. 23) we obtain 



q%)( V^-jf 0,0,0) <p(y)ii>($) + 



yean 



-£2 Mm 0, 0, 0) (V (y, 0)$(y, 0) + <p{y, 0)$'(y, 0) 
+ T (£2(^o) 1 ' + £ I (^/3o) r ' + 



+ 



(A.54) 



Using the expansion (B.116) of j'b and obvious relations for /3o(y; 0,0,0) and /3i (y;0,0,0), this 
expression is equal to 



q* (y) ( Vr-J^- P(P> v\v) v{y)^{v) 



yedn 



^ (3(0, y\y) (», 0)#(y, 0) + 0)$' (y, 0) 



4n 2 
ri/ / 1 + 6 
T I 3" 



+ 



/3(0,y|y)fc(g) + -J-^ /3(0,y|g) + 
n V 27Ti/ 



(A.55) 



+ C(r^) 



Using the normalization condition ( 2.1 5| ) concludes the proof of the expansion ( A. 46 ). By inspection 
we see that the /3-terms have form 



n V 2-kv 



2n 2 



6n 2 



(A.56) 
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B Geodesic theory 



Basic definitions 

In this appendix we review some facts from geodesic theory and list a number of useful expansion, 
some of which we have used in this work. The material related to a manifold without boundary is 
well known — see for example the classical works [l6| an d p7[. The theory of expansion near 
a boundary is less known. Some material can be found in ]18| |20|. The calculations are usually 
straightforward but cumbersome. We will present mostly only results. But all expressions presented 
were computed and checked by the author (see pi]]). 

We start with introducing the covariant expansion in a curved manifold. We would like to 
expand a sufficiently smooth tensor field A^^ on a manifold around a point x. First we change the 
dependence on a point z in the manifold M to the dependence on a vector Z from T x M, then we 
transform vector indices from different tangent spaces to one common tensor space, and finally we 
do the usual Taylor expansion of a linear-space-valued function on a vector space. 

To transform the tensor field on the manifold to a linear-space- valued function we need to know 
how to move tensors from one tangent point to another. We assume that we have given a metric g 
which define a parallel transport. It allows us to transform tensors from the tangent space at point 
z in a normal neighborhood of the point x to the space T x M along the geodetic joining these two 
points. In the normal neighborhood of x we can parametrize a geodesic by its tangent vector at x, 
i.e. we can define an exponential map u x 

u x :T x M^M , (B.l) 
ILR. Ux ( tX )=0 , ^-u x (tX)\ t=q = X . (B.2) 

If f(z) is some manifold dependent function, we use notation f{x;Z) = f(u x (Z)). This transfor- 
mation concludes our first step. Next we parallel transform vector indices of the tensor field to the 
space T x M along the geodesies starting from x. We define the tensor of geodesic transport t(x\z) 
from z to x 

b^ v (x\z) G T x M ® T* M , —l(x,u x (tX)) = , (B.3) 

or 

and its version with indices up and down 

I = i ■ g' 1 , L = g- 1 . (B.4) 

Using this tensor we can write down the tensor field A with transported indices explicitly. We obtain 
the linear-space- valued function on a linear space 

KZ{*\>) t a „(x\x;.)... i lv p(x\x;.)... : T x M -> T XJ * M . (B.5) 

Finally we can write the covariant expansion 

X;::-r:Z >,"„[;■ r: Z ) . . . l 1 v p(x\x\ Z) . . . = £ ^ . U^,, . . u- i /'< . . . /.>' . (B.6) 

feeNo 
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We call Ak(x) the coefficients of the covariant expansion of the field A at x. They are tensors at x 
symmetric in indices p u . . . , ii k . 

To compute these coefficients we need to develop geodesic theory to greater detail. First we 
define the world function cr(x\z) of the metric g. It is given by half of the squared geodesic distance 
between points x and z — see (1.18). For time-like separated points it is negative. The geodesic 
distance is then given by 

s{x\z) = \2o{x\z)\* . (B.7) 

We define geodesic tangent vectors er, a 

<j(x\z) — g~ l {x) ■ &\<j{x\z) , <j{x\z) = & r (j(x\z) ■ g~ l {z) . (B-8) 

Here, as before, di/ or d T f denote the gradient in the left or right argument of a bi- function f(x\z). 

The basic properties of the world function (see e.g. [^6|) are that its gradient vector a(x\z) is 
really tangent to the geodesic between x and z and it is normalized to the length of the geodesic. 
I.e. 

-Z =a(x\x;Z) . (B.9) 

We also introduce a special notation for the second derivatives of the world function 

a = VidiCT , a = did r <7 , a = V r d r cr . (B.10) 

To conclude our definition we introduce also determinants of t, i, a. They are well-defined objects 
— bi-densities on M 

i(x\z) = -^DeU(x|z) = Q~*(x)si*(z) , i(x\z) = 4r Det i(x\z) = q^{x)q^{z) , (B.ll) 

s(x\z) = \r>ct(- a (x\z)) . (B.12) 

Finally we define Van-Vleck Morette determinant 

A(x\z) = s(x\z)i 1 (x\z) = s(x\z) g^(x)Q^(z) . (B.13) 

For a bi-tensor F(x\z) on the manifold — a tensor object depending on two points in the manifold 
— we denote the coincidence limit 

[F](x) =F(x\x) . (B.14) 
The generalized Synge's theorem (see e.g. p^|) tells us that 

V[F] = [ViF] + [V r F] . (B.15) 

Coincidence limits and covariant expansions 

Equation ( |B.9| ) givesS 

a = ~ d z a ■ g- 1 ■ d r cr = i ditr • g' 1 ■ d x cr = ■ g ■ a = ■ g ■ a . (B.16) 

Taking repeatedly derivatives of this expression in both arguments we can derive 

a-g = a- a = a- a , g-a = a- a = a- a , (B.17) 
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and following identities 



cr^di M cr = 


2a 




(B.18) 


a Vi M cr = 
o- Vi M cr = 


* a 

a 


) 

; 


(B.19) 




G a/3 


- O a/Ji O p v g ^ -SXauPvV o , 




ff^Vi^cr^ = 
C M Vi M (7 a( 3 = 


O a/3 
C a/3 


- a aiJ.9' 1 **" & u/3 , 

— g~ * v a pa a v p , 


(B.20) 


CT^ViaVi^Vi^CT = 


C a/3 


- cr at j,ci p^g' 1 ^ u . 





Similar, more complicated relations hold for higher derivatives. Using the fact that coincidence limits 
of the world function and tangent geodesic vector are zero, taking coincidence limits of relations above 
and similar relations for higher derivatives, and using the Synge's theorem, we get 

M=0 , (B.21) 
[dicr] = [d r a] = , (B.22) 
[ViVitr] = -[ViV r o-] = [V r V r fj] = g , (B.23) 

[ViViViCr] = [VlVlVrff] = [ViV r V r Cr] = [V r V r V r cr] = , (B.24) 
[VlaVi/jVi^Vi^a] = -[Vi/jVi^Vi^VraCr] = [Vi^Vi^Vr^Vrca] = 

1 ( B -25) 

= — [VijyVr^Vr^Vro.cr] = [V r j,V r M V r pV r a cr] = — — (R aw a„ + Rai//3p) 

Similar relation for the fifth and sixth derivative can be found in pi] ] . 

Now we are prepared to compute at least some coefficients of covariant expansion. We start with 
the simplest case of the covariant expansion of a function / on the manifold. In this case we do 
not have problems with tensor nature of / and we do not have to worry about parallel transport of 
tensor indices. The equation ( |B.6| ) can be rewritten using (B.S) as 

E i - Z ^fk f , 1 ...^(x)a^(x,z)...a^(x,z) . (B.26) 
Taking derivatives of this equation and coincidence limits we can find that the coefficients are given 

by 

A,,., l =V ((ll ...V )li) / . (B.27) 

To do a similar calculation for a general tensor field A we need to know the coincidence limits of 
the geodesic transport tensor. They can be calculated from the equation 

a = -t-a (B.28) 

by taking derivatives and coincidence limits. We give only a list of some of them (sec |l7],[2^] or pi]]). 

\l]= 9 , (B.29) 
[Vii] = [V r t] = , (B.30) 

-[Vi/sViatpJ = [Vi Q V r /3i M J = -[V ra V r/ 3t M J = ^Ra/3/x^ , (B.31) 



2 

3 



[Vi-yVijOViaipi/] — — TjV-yRa/3/xiz — -V pRa-fni* , (B.32) 
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[ViaVi-yVi/jViat^J 



+ g R^5A RKa/3-y ff" lKA + g R/z^^A Rrect/3,5 .9 ~ 1kA + g R/^/3A Rkcx-^ 5~ 1kA + 

+ R/x^<5A Rre/3a- ) ' .9 + ^ Rp^-yA Rk/3q5 .9 + R/xi"5A .9 lKA + 

^ + g R^A/35 Rn/ta^ .9" K + g R^A/37 R 
+ g R^Act<5 Rk/x/3-)' .9" 1kA + g R^Act^ Rk^/36 .9 lKA + g R^Aa/3 Rr^iS .9 lKA + 
+ ^ R^re-y Rct/35A .9 lKA + ^ Rpii/K./3 Ra^SA .9 lKA + R/i^«/3 RaS-yA .9 lKA + 
+ ^ R/ii'KQ R/375A .9 lKA + ^ Rpii/rect R/367A .9 lKA + g Ri^/iKa R/3A7<5 .9 lKA 

Derivatives in other argument can be obtained using Synge's theorem and commuting covariant 
derivatives. 

Now it is straightforward to compute the coefficients in a covariant expansion of a general field. 
It can be done by taking covariant derivatives and coincidence limits of the rewritten equation (B.6) 

(-1)* 

feGN 



az:::(z)iZ(x\z)... c^{x\ z )... = £ L-L Ak %:: : ^ k (x)a^(x\z) ...*^(x\z) . (bm) 



We will not list explicit results. 

We can also expand a bi-tensor A(x|z) in both its arguments around some point y. We denote 
coefficients of such expansion A k j(y). I.e. 

t*(y\y; X)i*(y\y; Z)A(y;X\y; Z) = A k,l ^...^...u, (v) * Ml ■ ■ ■ ^ Z"* . . . Z v * , (B.35) 

fe,ieN 

where by L*{y\z)A(z) we mean a parallel transport of all indices from z to y. In the case of a bi-scalar 
f(x\z), similarly to ( B.27 ) we can derive that 

/mmi-Mh,"!-"! = lyp-Mi •••Vi Mj> )V(„/ 1 •••V rI/ ,)/] . (B.36) 

For calculations in appendix |X| we need the covariant expansion of the world function a(x\z). 
When we expand both its arguments at point y using the method described above we obtain 

a(y;X\y;Z) = ±(X-Zfg^(y)(X-Zy-±X»X»Z K Z x R l , KvX (y) + ... . (B.37) 

Clearly, the expansion of the world function at one of its arguments is given by equation ( B.16| ). 

Similarly, it is possible to derive (see |2l],|2^]) that coincidence limits of derivatives of the Van- 
Vleck Morette determinant are 

[A] = l , (B.38) 
[diA] = [d r A] = , (B.39) 

[Vx^Vi^A] = -[Vi M V rI/ A] = [V r/x V r „A] = -Ric^ . (B.40) 
and the covariant expansion 

A(y;X\y;Z) = l + ^(X-Z)-Ric-(X-Z)+... . (B.41) 
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Finally let us note that the Jacobian associated with a map 

m; 1 : z -» Z = -a(x\z) (B.42) 

is given by 

iDetDu^O)! = |Det(g _1 (ar) • a(x\z))\ = _1 (V) s(x\z) = i(x\z) A(x\z) . (B.43) 
(d — 1) + 1 splitting near a boundary 

Now we turn to investigate the domain ft with a boundary. We will study this situation locally 
— i.e. we will work on a neighborhood of the boundary with topology IxS where £ is part of 
boundary manifold. In such neighborhood we can perform a (d— 1) + 1 splitting which is discussed 
for example in |23| . It is given by a time function t and time How vector t such that t ■ dt — 1 . We 
use notation of usual 3 + 1 splitting of spacetime even if we do not necessarily assume that t plays 
the role of a time coordinate. We assume that the condition t — defines the boundary and that 
t > inside of the domain f2. We denote E t hypersurfaces defined by conditions t = const. We 
denote n and n inside oriented normalized normal form and vector, q orthogonal projection of the 
metric g on the hypersurfaces E t , and b orthogonal projector to hypersurfaces St. I.e. 

g = n 2 nn + q , g' 1 = n nn + q' 1 , S = iin + b , (B.44) 

where q' 1 is inverse of g in the tangent space to the hypersurfaces. The phase factor n governs the 
signature of the metric g and the character of the hypersurfaces. We will use shorthands 

^:::'r.v x - = r? n IL . (b.45) 

We also use 

-* dcf -» <-» def <-> «- del «- /„ ,„n 

= GTii II , = o-,m, , cr„ = cr M „ . (B.46) 

Decomposition of the time flow vector t defines lapse N and shift N 

t = Nii+N , dt = Nn . (B.47) 

We denote d the hypersurface gradient — an orthogonal projection of a spacetime gradient to 
the hypersurfaces E t 

D/ = b-d/ , (B.48) 

and V the hypersurface connection of the metric q. It is related to the spacetime connection as 

VA = b*VA for A such that A = b*A . (B.49) 

where by b*A we mean projection of all tensor indices to the spaces tangent to the boundary. We 
denote by R, RIC, TZ, and V 2 the Riemann curvature tensor, Ricci tensor, scalar curvature and 
Laplace operator of the metric q. 

The extrinsic curvature K is given by covariant derivative of the normal form 

K = b ■ Vn , (B.50) 

and we use shorthands 

k = K uv g 1 ^ , 

, , (B.51) 
K 2 = K Kfi K Xv g lKX g- 1 ^ . 
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We define the time derivative of a tensor field A tangent to the hypersurfaces: 

A' = b*C r A for A such that A = b*A . (B.52) 

Now we list a number of useful relations between spacetime quantities and "space" quantities, 
derivations of which are straightforward and for the case n 2 = — 1 can be mostly found, for example, 
in [El. 



n ■ Vn = — D In N , 
V • n = n 2 k , 

n-Vq = n(n (Din AT) + (d In N) nj 

&7 V M g a/3 = -n a Kp^ - npKay , 
V • b = -n 2 nk + D In N , 
b • (Vb) • b = -Kb 
q = 2n 2 NK + £$q , 



(B.53) 
(B.54) 
(B.55) 



K' = TV (iJ • VK) |m| + 2n 2 NK ■ q' 1 - K + C^K , (B.56) 
k' = NS-dk + Nk . 

The curvature tensors of the spacetime metric g and of the space metric q are related by 

Rnct n/3 117 ii(5 = Ra/3-y(5 + n 2 (Kq^K^ — Kq^K^) , 
R|iau/3T_L = ™ 2 (V Q K(3 7 — V(3K a7 ) , 

R ll±ia = n 4 (Kfc-K-g 1 • K) - n 2 ( VV In AT + (d In AT) (d In N) + (V • (HK)) MM ) = ( ^ B " 57 ' ) 

= n^K-q- 1 -K-n 2 (vv\nN + (d In N) (Din AT) + — K' - — £^K) , 
Ric„„ = Ric - Win AT- (Din AT) (Din AT) - (V • (HK)) , 
Ric, a = n 2 (v ■ q' 1 ■ K - uk) , 

Ric ±x = n\k 2 -K 2 ) - n 2 V • (nk + q 1 • (Din AT)) = (R58) 
= -n 4 K 2 - n 2 ( V 2 In AT + (d In AT) ■ q" 1 • (d In N)+n- dk\ 
R = Tl + n 2 {k 2 - K 2 ) - 2V • (iik + q- 1 ■ (Din A^)) 

= 11- n 2 (k 2 +K 2 ) - 2 f 12 • dfc + V 2 In A^ + (d In AT) • q 1 • (Din AT)' 



(B.59) 



Geodesic theory near a boundary 

We can develop geodesic theory on a hypersurface £ similarly to what we did for the spacetime 
M. On the boundary we denote the exponential map v x , the tensor of geodesic transform and 

its determinant q and j, the world function, its derivatives and its determinant p, p, p, p, p, p 
and r, and Van-Vleck Morette determinant A. Finally, we denote {.} the coincidence limit on the 
boundary. 

In the neighborhood of a part of the boundary £ of the domain f2 in which geodesies orthogonal 
to the boundary do not cross we can also define the map w 

w:Exl->M , 

(B.60) 

w(x,rj) is geodesic , w(x,0) = x , w (x,0) = n . 
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It maps point x on the boundary "orthogonally" to the domain 17 by the distance r\. We denote E^ 
the hypersurface which we obtain by shifting E = Eo by the distance r\. We also use the notation 

w n : E — > , w v (x) = w(x, r/) , 

_ x (B.61) 

This foliation is special case of the foliation discussed above. We obtain it for the choice of lapse 
and shift N = 1 and N = 0. 

For a tensor field A(x) on spacetime we denote by A(x, its dependence on x and £, and ^4(2, 
the tensor field on the boundary manifold obtained by transformation of tensor indices of A(x,£) to 
the boundary tangent bundle using u>_^ 

A(x,0 = A(w(x,0) , A(2,0 =t&V(^(2)) , (B.62) 

where is the induced transformation on tangent bundles. For a bi-tensor A(x|z) we mean by the 
boundary coincidence limit 

{A}(y;Z,Q=A(y,Z\y,0 ■ (B.63) 

Specially, we have a metric q{y,vi) (generally different from q(yj) on the boundary manifold, 
volume element q^(y,rj), and associated connection V. It is related to the connection V by 

V = V © 7 . (B.64) 
The relation of corresponding curvature tensors is (see e.g. p3|) 



R 7C //3 = R 7 c//3 + V 7 7i/3 - + it A - 7^7^ 

KlC a p = KlC af} + V M 7^ - V a7 ^ + 7a M 7^ " 



(B.65) 



From the definition of the map w we have 



{^K9;£,C) = (£-0 2(5,0 . (B.66) 

{d iC t} = . (B.67) 



(B.68) 



Differentiating this equation we obtain the differential map Dw 

Bw u (x) :T x E£->T a E c , z = & 5 , c (:e) , 
D^w iX (x) = - (Vi M Vi K (r) (x\z) a~^ KV (x\z) . 

In the special case £ = we get 

^w v {y) = -( W i 2 K^(y) + a^ n {m,l)) a^ KV {y\y,rf) . (B.69) 

Here c' 1 is the inverse of er M in spaces tangent to hypersurfaces £,,. Because Dto^ = Dii^ ^ we 
have 

{^,^} = {(vM( v -V r Aa)^ lKA } . (B.70) 
Using this relation, the definition of the Van-Vleck Morette determinant and 

n(x) ■ a(x\z) = — n(z) ■ g(z) for z = w^^(x) , (B.71) 
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we get an expression for the Jacobian associated with the map w^^, 

j(y;£,C) - |Det„D© € , c | = 

= {A" 1 q" 1 (Det„ Vi V 1( r) } (y; £, Q = {A q(Det„ V r Wr)" 1 } (y; £, C) 

As special cases we have 

This also gives the expression for the Van-Vleck Morette determinant, 

i 

A(y, elf/, = q 4 (S, £) (y, ((Det„ ViVkt) (Det„ V r V r fi)) \y, £\y, C) 
Finally we can prove that 

= q\y) (Det„ ViViff)(fJ, CIS, = qXv) (Det„ V r V r ?)(y, C) 



(B.72) 



(B.73) 



(B.74) 



(B.75) 



(B.76) 



Thanks to flB.67|) , (|B.72|) , and QB.73D , the function j defined in ( |A.36| ) is 

Km L = j(y, A(g, £\y, 3(y, ^ (y) (Det,, v r v r( r) 2 (g, c) = 

= j(y,OA(y,C|j/,C) q%0 (Det„ V r V rt r) 2 (y,£|y,C) = 

= j(as)A(x|a!) q*(«) (Det„ V r V r a)^(x\z) = j(x) A*(x\z) J(y;^O h 

= (](y,0kv,m,0j(y,0) i ■ 

The equation previous to this one is a straightforward consequence. 



Reflection on the boundary 

In section we have worked with the geodesic reflected on the boundary. We recall its definition 
here and list some useful properties which allows us to prove the relation (A. 52). 

We will study the geodesic Xb{x\z) between points x and z which is reflecting on the boundary at 
point b(x\z) — an extreme trajectory of the functional given by the half of squared length, with the 
condition that it has to touch the boundary. We use the convention, that for any quantity depending 
on two spacetime points f(x\z) we denote 



Mx\z) = f(x\b(x\z)) , Mx\z) = f(b(x\z)\z) 



(B.77) 



If we denote the parameter at which the geodesic reflects on the boundary A r (x|z) and its complement 
Xi(x\z) 



b(x\z) = x h (x\z)\ Xl ( x \ z ) € dil 
1 = Xi(x\z) +\ T (x\z) . 

we can write the reflected geodesic as joining of two geodesies 

[r,Xb] = [Ait,*i] [A r r,S r ] 



(B.78) 



(B.79) 
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The extremum conditions on position of the reflection point and reflection parameter are 



(do-)i (Dcr) r 



Ai 



^L = ^L 



where D acts in the argument on the boundary. 
We define reflection world function at, as 



Clearly 



\\ A r Aj^ Aj? 



Ai — 4/ — > A r — a I — , 



CTb 



CTb 



Using the last equation we can get 

dlV^b = (diVff)l , d r y/<Tb = (drVcr)r 

Similarly to the case without boundary we define 



Sh 



^(Tb = Si + S T 



-» -1 , Cl <- _i 

<?b = g -diC7b = -r" ' cr b = d r Cb-g = -— , 
A] A r 

Cb = ViViCTb , o"b = did r <Tb , crb=V r V r (7b ■ 

Additionally we define 

1 1 
s± = -— n(b) ■ (d r er)i = -— 5(6) • (dicr) r = n ■ (di cr b ) = 3 • (d r o- b ) 
M A r 

and we denote differentials of maps a; — > b(x\z) and z — » b(x\z) as 



(B.80) 



(B.81) 



(B.82) 

(B.83) 

(B.84) 
(B.85) 
(B.86) 

(B.87) 
(B.88) 



i.e. Jf we displace points x and z in directions X and Z, the reflection point moves in direction 
X ■ b (x\z) + b (x\z) ■ Z. Finally we define the reflection Van-Vleck Morette determinant Ab 



A b = |Det CTbl i" 1 



(B.89) 



Some long algebra gives 

Ob 



- b ■ B ■ b - — !— ^(dicr b ) (di cr b ) + , 
Za b Ai Ai 

<-> — 1 

cr b = - & • B ■ b + - — (dicr b )(d r (Tb) , 
Za b 

a b = - b ■ B ■ b - — !— ^-(d r (Tb)(dr(Tb) + , 
ztTb A r A r 



(B.90) 



where 



B = 2^ (v r V r V^)i + ViViv/o) r ) 



<T M 1 (T || f 1 (Tb 



Ai 



2 (Tier, 



(D 1 a) r (D r a) 1 + 2s ± K(6) . (B.91) 
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Using these relations, a more intricate calculation gives the space coincidence limits 

{£ bll } = -{ 6 • B ■ X } , 

, {%}= , 

{(ViVidb)" 1 } = {a h \ ■ (VrVrCTb)" 1 • CT~J,} = 

= ■ (V r V r a)i • (A^VrVrCr^ 1 + A r (Vi Vict); 1 ) • (V r V r cr)i • a-\j , 
{(ViViO-b)" 1 • a bn } = {a]l • (V r V r cr b )} . 



(B.92) 
(B.93) 

(B.94) 
(B.95) 



Here inverses are taken in the spaces tangent to the boundary. Taking the determinant of the last 
equation, we find 



{(Det„ CT b „) 2 } - {(Det„ ViV x <r b ) (Det„ V r V r <r b )} , 
A b (y, = f % q" 1 % (Det„ Vi Vi<r b ) (Det„ V r V r a b )) 2 (y, 



(B.96) 



and 



7 2 fx) 7 2 fz) 

J —^- (Det„ ViVitrb)(a;|2) - ^7-^ (Det„ V r V r o- b )(:c|z) 
qw qO) 



q(y) 



(AiCVpVptr)! 1 + Ar(ViVio-); 1 ) (x\z) 



(B.97) 



for a; = w(y, £) and z = ■&(?/, £). Putting these relation together we obtain 



j(x) A h {x\z)j(z) 



q(x) 



(Det„ ViVi£Tb)(a 



q(*) 



(Det,, V r V r cr b )(a;|2;) 



(B.98) 



for a; = -&(?/,£) and z = w(y, C)- The equality in ( A. 52 ) is a straightforward consequence of this 
relation. 



Covariant expansions near boundary 



Finally we will write down coefficients in covariant expansions (A. 32) and (A. 33) of the world 
function a and function I defined in (A.3C). These are expansions inside of the boundary manifold 

of to-mapped functions cr(i,£|2,£) and l(x,£\z,Q around point 2. The derivation is long and 
technical. It uses the general method discussed above and a transformation of the connection V to 
the connection V. Fortunately, we need only the spacetime coincidence limit of the coefficients (i.e. 
&k,i{y\ *7, which simplifies the calculations significantly. But even then the calculations is too 
long and uninteresting to be included it here. We list only the results. See also [p~8H20|l for similar 
calculations. 

The coefficients of the boundary covaria nt exp ansion of the spacetime world function a at some 
general point y (slight generalization of eq. ( A. 32] )) are 



fro,o(9;£,C) = 2 n2 (£-0 2 . 
&o,i(i/;£>C) = 3i,o(y;£,0 = , 
[£2,0] = -[STi] = [^0,2] = q , 

^(a(3^-f)fJ. ' [^2, lot/3*] = [?l,2(ta/3] 



~ /x ~ 



(B.99) 

(B.100) 
(B.101) 
(B.102) 
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[<?4,0 a.f3f6\ — [9o,4a/37i] — 
[&3,1 af3-f k] = [^1,3 k a/37] = 

[52,2ci/3ka] = — - (Rctft/3A + R«A/3k) — 



(B.103) 



1 



1 



-n z K aK K px - -tj'KcaK^ - n K a/3 K KA - 7a/37 K A§/«' • 

The coefficients of the b ounda ry covariant expansion of the function I at some general point y 
(slight generalization of eq. ( A.33f )) are 



[?o,o] = 21nj , 

[h,0a] = [1 0,1 a] = lay, ' 
['2,00/3] = [J 0,20/3] = 

= \ (-Kp + 2n 2 K a , l Kf, u q~ 1 - n 2 K af3 k 



fl 1,1 al 3} = l(K 



Computing normal derivatives we also get 

1 



n 
1 



(In j)' = nk , 
r(ln])" = fc' , 



and 



-[(i^-) r ']-o , 

^[(lnj) r ' r ']--J(fc'+n 2 K 2 ) . 
Finally, we have boundary coincidence limits 

{<M(y;£,C) = ^ 2 (£ + 2 , K}(S/;f,C) = ™ 2 (£ + , 
{2 b }(B;&C) = te + CMfcO , {^b}(y;^C) = K + CKy.O 



{Ai}(y;£,C) = 



> {A,}(9;e,0 



Using these relations, equations (B.92), (B.93), with help of (B.6E) and 



we can derive 



-{vbnW, £, C) = q(y) + (( + (- 2-^-)n 2 K(y) + + C) 2 ) 



(B.104) 
(B.105) 



(B.106) 



(B.107) 



(B.108) 



(B.109) 
(B.110) 
(B.lll) 

(B.112) 
(B.113) 
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From this follows 



Together with 



{A,.}(//:i.O I : (e + C-2^)n 2 fc(y) + O((e + 2 ) 



(B.114) 



](y,r]) = l + n 2 k{y)r ] + 0(T 1 2 ) 



(B.115) 



and the definition (A. 52), it finally gives the expansion for j'b, 



Jh (y; t,\,Q = l+(i±i-(l- 2p)X ) n'k(y) + 0((£ + ( + Xf) 



(B.116) 
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C Special functions 

In appendix [X] we have used various integrals of exponentials of quadratic exponent and integrals 
of such integrals. Here we summarize properties of such integrals. We will introduce a special 
function R y closely related to error functions erfc(x), the definition and properties of which can be 
found, for example, in [p4f . The derivation of the properties below are not all simple, and we do not 
include them here. 

We define the function K u (x) for positive v as 

R v (z) = f ^-r J dxx^exp^x-z) 2 ) . (C.l) 

It is a solution of the differential equation 

R' u (z) = vR v+1 (z) - zRv(z) , R„ *~ > ~°°> . (C.2) 
In the limit v — > and for v — 1 we have 

RoO) = expf-iz 2 ) , (C.3) 

(C.4) 



R 1 (^) = V2^-y|erfc(-^) 



We also have the recurrence relation 

Ru+2{z) = -—j(zR 1/+1 (z) + R v (z)) . (C.5) 

For ceN these functions are combinations of Ro and Ri with polynomial coefficients 

R«+i = P«Ri + q«-iRo for neN , (C.6) 

where 



Pn+l(*0 = —(zp n (z) + p n -i(z)) , po = 1 , Pi = Z 

n + 1 

q n +i(z) = — —(z<\ n {z) + q n -i{z)) , q = 1 , qi = ^z 



(C.7) 



These polynomials satisfy 

p', = Pn-i , q' t = (n + 2) q„+i - p„+i , (C.8) 

and 

V27p„(z) = R n+1 (z) + (-l)"R„ +1 (-z) = ^^-H„(^) , (C.9) 
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where H n are the Hermite polynomials (see p4fl). 
Values at zero are 

a r(§ + 1) 1 V2tt J ( 1/ _ 1 1 )M f° r v natural and even 

^ 2 r(z/+ 1) ~ 22T^r(^- + 1) ~~ 1 for i/ natural and odd 



J -TT for n even , ) for n even 

Pn(0) = < P„ = < ! > 

I for n odd I ^ n _ 1 ^ ! , tor n odd 

q B (0) = (^ I J ir forneven q'={° for n even 

\° for n odd ' " - (w+ 1 1)!! for n odd 



The behavior for small z can be found for natural v with help of relations ( |C.5| ) and 

Row = E t^ 2fc 



(2fc)!! 



V fcGN 



( 1 ) * - z 2 ' £ + 1 



(2fc+ 1) (2ft)!! 

The behavior for \z\ ^S> 1 and n £ No is 

R n+1 (z) = V2^p„(z)e(z)+exp(-iz 2 ) O(^pr) , 

where is the step function. 

Now we can write down results of some integrals in terms of these functions. For n £ N 



— / dxx n exp( (x — z) 2 \ = s/Ztt p n ( z ) 

n\ J V 2 / 



Further for ft, I £ No we have 

\l+m 



/ \ m—0,...,l 

{ — oo,x) 



k+l-2m 



i— t w , x \/2^ (2m-l)!! //c-l N 

2m^fe-l 

E -V n iW™ +1 (-z) 



i 

{0,x 



(-l) fe R fc+(+ i(0) + yr^-TTT E 



m! 
m=0,...,i 



2^ ^ (2m-l)H/'fc-l^ a;fc+I _ am 



H(fc-l)! ^ k + l-2m\ 2m 

2m<fe-l 



E A« m R* + ,-^ 1 (-a:) 
* — » m 



m 

m=0,...,Z 



Finally for to, ft, Z G M + and n = m + ft + I we have 



1 ' dt d C + Cf) = ^ 2-^ F(m + * + 1) IV + Z + 1) 



n! 7 s s (£ + C) m V 2 VS S W r(6 + m + l)r(2^±i) 
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Notes 



1 We use the MTW sign convention (p3|) for spacetime quantities generalized to accomodate 
both Euclidian and Lorentzian versions of the theory. The version is determined by the con- 
stant factor n which is equal to 1 in the Euclidian version or to i in the Lorentzian version. 
Specifically, geometric quantities defined in spacetime M with an embeded hypersurface E are 
related in the following way: 

metric: g a p = n 2 n a np + q a p , q positive definite , 

volume element: = — (Det g)? , 

curvature: Rap^a" = V[ a V /3] a' J ' , Ric Q/3 = , R = Kic af3 g~ laP , 

extrinsic curvature: K a p = b^b^ V^n^ , k = K a pg~ laf3 , 

n has the orientation of the hypersurface £ , 
projector on £: = & n a , n a = n 2 n^ 1 ^ . 



2 Euclidian and physical actions are related by 

-7(h) = uS(h) . 

The version of the theory is determined by actual value of the factor v — whether it is real or 
imaginary. 

3 The factor n, which governs the signature of the spacetime metric, is choseruhere for conve- 
nience and reflects that we are using Lorentzian convention for volume element!! I.e., the phys- 
ical amplitude is —K, but the quantity K will have nicer properties in language of Lorentzian 
quantities. Similarly, for the Green function the physical amplitude is —G F , but we will use 
often the quantity G F to express properties of the amplitude. 

4 To abbreviate formulas we often use different dots to indicate contraction in different vector 
spaces. We use "•" for contractions of tangent tensor indices, i.e. g a p aP — g ■ a. We use 
"•" for contraction in functional vector spaces of function and densities on spacetime, i.e. the 
bullet means an integration over spacetime (or an action of a distribution on a test function 
— a formal integration): 




Similarly we use "■" for integration over a hypersurface E or a boundary dfl of a spacetime 
domain. 
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5 A volume element on spacetime defines a bi-distribution Q = — a delta function 
normalized to the volume element, i.e. <j> • Q • ip = J (j>tp g5 . Clearly for a smooth function / 
we can define a distribution fQ. 

It is convenient to represent differential operators on the manifold M as bi-distributions. We 
use arrows <~ and ~> to indicate direction of derivatives. So, for example, 

tp • (d a a a ) »u = uj • (a a d a ) • ip = J uj a a d a tp 

for a test function tp, a test density uj, and a vector field a. 

6 A meaning of the tilde in the index of the quadratic form Ty. will not be discussed in this 
paper. It will be explained in the following paper and can be safely ignored here. 

7 Again, we factorize out the prefactor n motivated by the fact that J is a density, i.e. propor- 
tional to volume element jp. 
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